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Abstract 

The transverse momentum spectra of secondary 77 particles produced in P + P, D + Au 
and Au + Au interactions at y^sjvTV = 200 GeV at different centralities have been studied 
in the light of a non-extensive thermodynamical approach. The results and the possible 
thermodynamical insights, thus obtained, about the hadronizing process have also been 
discussed in detail. 
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1 Introduction 



The multiparticle production in relativistic heavy ion collision provides important information 
about the hadronizing process in the aftermath of any high energy collision phenomenon. Vari- 
ous phenomenological/mathematical models have, so far, been utilised to analyse the available 
experimental data to understand the different characteristics and dynamical aspects of this 
process. Amongst the various existing models, applicability of Tsallis non-extensive thermo- 
dynamical approach [TJ [2} [3J HI EJ El [7j , which has already left it's marked impression in other 
branches of physics, has recently gained pace in understanding mainly the statistical as well 
as the thermodynamical characteristics of the production of various secondaries in high energy 
nuclear interactions at different energies S[ini[ni[I2[]2[]Il[]i2 

EU ESI EHl E3 EH]- We have taken up, here in this paper, the task of analysing the transverse 
momentum spectra of one of the heavier secondaries, //-mesons, produced in different nuclear 
reactions at RHIC-BNL in the light of Tsallis non-extensive thermodynamics. This is in con- 
tinuation of our previous work[29j, wherein transverse momentum spectra for the secondaries 
like 7r + , 7r~, P and P were treated with the same approach. As in our previous work[2S], here 
also we have confined ourselves only to a particular interaction energy {s/snn — 200 GeV), so 
that a comparative study can be made among various nucleus-nucleus systems. Our mention 
of the P + P-interactions in the same breath and bracket with the host of nuclear reactions 
is to be viewed somewhat loosely; this was done with our eyes cast on the same interaction 
energy alone. This is specially so when the 'centrality '-aspects are touched upon and dealt 
with. In our actual work-procedure we have been attentive to the particularities of each in- 
teraction on a case-to-case basis. In the recent past, Biro et al |19|, [20] analyzed successfully 
transverse momentum spectra of different identified hadrons, including ?7-mesons, produced at 
RHIC energies, on the basis of both the non-extensive statistical approach and the very essence 
of quark-coalescence model. However, the study was, probably by choice, confined to the min- 
imum bias Au + Au collisions in central rapidity region. Our present study, however, takes 
into account (i) the effect of entire longitudinal momentum^) range, (ii) the role of the same 
in different nuclear collisions(Aw + Au, D + Au and P + P) and (iii) the impact of centrality 
dependent windows. 

The r] mesons are the heavier species among the various mesonic secondaries produced 
in high energy nuclear reactions. The formation of such heavy partonic bound states and 
their interactions in nuclear matter under different conditions throws light on the nature of 
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the possible hadronic phase transition to a partonic state, with a greater degrees of freedom, 
formed immediate after the nuclear collisions at relativistic energies. The transverse momentum 
spectra of such secondaries at different centralities provides significant clues on the mechanism 
of the formation of such bound states, on their partonic constituents as well as on the thermal 
and dynamical characteristics of the regions in the partonic matter where they emerged from. 

The work is organized as follows: Section 2 presents an outline of the hadronization process. 
In Section 3 we give some very brief sketch of the nonextensive statistics and the main working 
formula to be used in our work here. The results are reported in next section(Section 4) with 
some specific observations made. And the last section is preserved for our conclusions. 

2 Non-extensivity and Hadronization Process: An Out- 
line 

The evolution of the partonic system created in RHIC-BNL experiments is generally believed 
by a large section of the theoretical physicists to be best described by hydrodynamics of an 
"almost ideal fluid". This approach really gives a fair description of data on the particle 
transverse momentum spectra and asymmetry of the transverse flow. However, the prospect 
of this theory and methodology is, at the same time, riddled with and marred by two very 
serious problems: (i) the 'HBT puzzle' and (ii) the 'puzzle of early thermalization'[30| [31] 
which we are not going to elaborate here in detail. Besides this, the proposed microscopic 
mechanisms lack in the detailed knowledge of fragmentation and recombination etc. These 
deficiencies have prompted the physicists to turn to many different models, often from various 
different angles, of which the non-extensive thermodynamic approach is an effective one which 
is, thus, essentially a thermal model. This model assumes the formation of a system which is 
thermally and chemically in a near-equilibrium state in the hadronic phase and is characterised, 
in the main, by two thermodynamic variables(observables), for the hadronic phase. These two 
are 'temperature' and 'non-extensive entropy' of the system attained finally. The deconfined 
period of the time evolution dominated by the constituent partons(quarks and gluons) remains 
hidden: full equilibration generally washes out and destroys large amount of information about 
the early deconfined phase. Successes of thermal statistical models are shaped and conditioned 
by the blackening out of information on some intermediate states. And this is an accepted 
reality. 
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In the case of full thermal and chemical equilibrium relativistic statistical distributions are 
generally based on exponential nature of the spectra for the pr- distribution of the produced 
hadronic secondaries. But, observationally experimental data at SPS and RHIC energies reveal 
strongly the traits of non-exponential and power-law nature at high-p^ which is also supported 
by the dictates of the perturbative QCD[32, [33], |3l]. So, for the overlap or the transitional 
regions of the pr- v alues(low-to-high), a stationary distribution of strongly interacting hadron 
gas or quark-matter in a finite volume could be considered to have such a distribution with dual 
nature. In fact, Tsallis distributions 3, which we have resorted to, satisfy such a criterion[36] . 



3 Nonextensive Statistics and Transverse Momentum Spec- 
tra 

The nonextensive entropy by Tsallis generalized statistics is given bypQ, 

s q = ( i - £ pD (!) 

where are probabilities associated with the microstates of a physical system with nor- 
malization YU Pi — 1 an d Q is the nonextensivity parameter. For q — ► 1, eqn.(l) gives the 
ordinary Boltzmann-Gibbs entropy 

S = - £ Pt ln Pl (2) 

i 

The generalized statistics of Tsallis is not only applicable to an equilibrium system, but 
also to nonequilibrium systems with stationary states[9j. As the name 'nonextensive' implies, 
these entropies are not additive for independent systems. For a system of N independent 
particles, where particle 1 is in energy state e^, particle 2 in energy state e i2 , and so on, the 
Hamiltonian of the system in nonextensive approach is given by [9], 

1 However, the coinage 'Tsallis distribution' is not fully justified, because Vilfredo Pareto[55] used much 
earlier and long ago a power-law probability distribution found in a large number of real-world situations, 
especially in the field of Economics. And beyond the borders of economics, it is generally referred to as the 
Bradford distribution, though the high energy physicists, since 1970s, prefer to term this as a "cut-power- 
law distributions" . In the context of non-extensive thermodynamics alone, this is normally familiar as Tsallis 
distribution. 
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H(i 1 ,i 2 ,...,i N ) = e h + (<?-!)/? Yj % e h + - l) 2 /? 2 Y UiU&t + ■■■ (3) 

3 j,k j,k,l 

where (3 = 1/T is the inverse temperature variable. The above equation clearly indicates 
that for a non extensive system the total energy is not the sum of the single-particle energies. 

The energy associated with a particle, which is a hadron in the present context, denoted 
by j in a momentum state % in a fireball produced in a high energy nuclear collision is given by, 



e<j = \jv1 + rrij (4) 
The nonextensive Boltzmann factor is defined as [9] 



xn = (1 + {q-l)^- qKq - X) (5) 



with q — > 1, the above equation approaches the ordinary Boltzmann factor e - ^. If z/, 
denotes the number of particles of type j in momentum state i, the generalized grand canonical 
partition function is given by, 

z = EIK* ( 6 ) 

The average occupation number of a particle of species j in the momentum state i can be 
written as[9j 

Vij = -'''^ = (1 + (q ~ ± 1 

where — sign is for bosons and the + sign is for fermions. 
The probability of observation of a particle of mass mo in a certain momentum state can be ob- 
tained by multiplying the average occupation number with the available volume in momentum 
space j9]. The infinitesimal volume in momentum space is given by 

dp x dp y dp z = dpL pt sinO dpx d6 (8) 



where pr = yp% + p\ is the transverse momentum and p z = pl is the longitudinal one. 
Hence, the probability density w(pr) of transverse momenta is obtained by integrating over all 
6 and p^\ 
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/+oo \ 
dph , (9) 
(1 + (q - 1)13 J p 2 T + p\ + ml)i/^- 1 ) ± 1 



Since, the integrand in the right hand side of the above equation is an even function of pi, one 
can easily write, 

r+oo \ 

w(j>t) = const, px I dpL (10) 

JO h ±(n - 1 \R< /T)2 _|_ 4- m.2 W(g-l) 4- 1 



1 + ( q - l)^p 2 T + p 2 L + m^/t'- 1 ) ± 1 
Since the temperature T, which will be called here as Hagedorn temperature and will be 



denoted as T e ff, is quite small, i.e., /3ypj, + p 2 L + ml 3> 1, one can neglect ±1 in the denomi- 
nator of the previous equation. Hence, the differential cross section, which is proportional to 
w(pt), can be expressed as, 

1 1 (iff f°° , . (q 

oc " ■ ' ' ' 



a px dpT 



rdp L (l + ^ml + pl + pl) (11) 

Jo T eff 



or, in terms of the transverse mass(m? r = wip + pj,) of the detected secondary, 
1 1 da f°° (q 



. oc / dp L (l + ^—±^m 2 T +pl)-^-V (12) 
amTdm T Jo T e ff v 

If we put x = -P^- and u = ^ aL , the above equation can be written as, 

^ Teff T ef f ' ^ ' 



1 da r°° , , , / x 2 

oc u 



POO / rpA 

A &(l + (g-lWl + -)-^- 1 ' (13) 

a drriT Jo V u 

For large x the integrand is small and it's contribution to the integration can be ignored[8]. 



Hence, for large u and small x we can assume, + ^ ~ 1 + J^-. The inclusion of this 
assumption in the above equation yields 

' ^ oc u(l + (q- l)u)-^-V r dx( 1 + - / ( ; q ~ 1)X \ ; (14) 



a drriT Jo 2u(l + (q — l)u) 

This form of equation is in accordance with the right hand side of eqn(28) given in [8]. So, one 
can write, following the eqn(32) of the same reference, 

I_^L oc u 3 / 2 (l + (q-l)u)-^ +1 * (15) 
a drriT 

or, 

-4^~ = ex w 3/2 (l + {q- l)u)-^ (16) 
6 



where c\ is a normalization constant. 

The average multiplicity of the detected secondary in the given rapidity region can be 
obtained by the relationship 



< N >= - / dm T = d / u 3/2 (l + (q - l)uy— + 2dm T (17) 

<T J mo Ctni r p J mo 

Hence, the constant c\ can be expressed in terms of < N > by the relationship 

<N> 

Cl = a 1 18 

%c« 3/2 (i + (?-i)«r" + ^« 

where u = -S^. 

Combining eqn(14) and eqn(16) one can write 

1 da <N> o /2/i . , 2_,i 

= „ — u 3/2 {l + {q-l)u) 19 

a dm T Teff j™ u 3/2(! + ( 5 _ +2rf M 

Since, 

1 1 da 1 1 da dN 

a m,T dm T a p T dpT Pt dp t 

u 1/2 (l + (q-l)u)-^i + -z (21) 



we can write 

dN < N > 



Pt dp t T 2 ff u 3/2^ + ( ? _ i) M )-^+2rf n 
Equation(21) provides the working formula for the present analysis. 



4 Results 

The working formula was applied, in it's present form, to obtain a fit to the data on r\ production 
in P + P collision[Fig.l(a)] at ^/snn = 200 GeV and the corresponding parameter-values are 
given in Table-1, where n Q denotes the average multiplicity of r\ produced in P + P interaction. 

However, when the same formula was employed to analyse data from nucleus-nucleus 
collisions like D + Au and Au + Au collisions, the parameters, which had been set free for 
P + P collisions, were constrained by the following relationships [24J: 

T eff = T (l-c(q-l)) (22) 
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< N > -n N part 
<N> 



c(g-l) 



(23) 



with c 



where D, c p , p, T are respectively the strength of the temperature fluc- 



Dc p pT 



tuations, the specific heat under constant pressure, density, the temperature of the hadronizing 
system when it is in thermal equilibrium(g = 1) and N part is the number of participant nucle- 
ons. Eqn.(22) describes the fluctuation in temperature where it is assumed that the effective 
temperature T e ff is the outcome of two simultaneous processes: (i) the fluctuation of the tem- 
perature around T due to a stochastic process in any selected region of the system and (ii) 
some energy transfer between the selected region and the rest of the system, denoted by 0[24|. 
It is absolutely uncertain whether such energy-transfers could/should be invariably linked up 
with flow-velocity(normally denoted in the hydrodynamical model-texts as 'it'). So, for the 
sake of calculational simplicity and correctness we assume the factor 0, for the present, to 
be independent of any flow- velocity. The fluctuation in multiplicity is described by eqn.(23). 
The assumption behind this relationship is that if N-particles are distributed in energy accord- 
ing to Tsallis non-extensive distribution, then their multiplicity will obey Negative- Binomial 
distribution [24|. 

The fits for D + Au and Au + Au collisions obtained on the basis of equation(21) alongwith 
the constraints given in eqn(22)-(23) are depicted in Fig. l(b)-Fig. 1(d). The values of various 
parameters obtained from the fits are given in Table-2. The values of T e ff and q calculated from 
the fitted parameters are given in tabular form in Table-3 and in graphical format in Fig. 2 as a 
function of participant nucleons. Fig. 3 depicts graphically the behaviour of average multiplicity 
of rpmesons produced per pair of participant nucleons as a function of N part . The normalized 
values of < iV > for different nucleus-nucleus collisions at ^/sNN = 200 GeV remain almost 
constant with respect to N part which indicates the linear dependence of < iV > on N part , and 
hence on the system size. 

The obtained values of q and T e ff, on the average, show centrality-dependences for almost 
all the collisions. The values of To and c obtained from different fits show almost constant 
behaviour according to expectation and the average values of them are found to be 155 MeV 
and 1.88 respectively. The performance of the present non-extensive approach vis-a-vis the 
studied experimental data could be rated to be moderately satisfactory, as is seen from the 
X 2 /ndf- values given in the last column of Table- 1. 
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5 Discussion and Conclusions 



The centrality dependences of the parameters q and T e ff are quite clear from the Fig.2. The 
systematic trend of the non-extensive parameter q, in general, is that it decreases with increase 
in the number of participant nucleons, i.e. with the increase in centrality. On the other hand, 
the Hagedorn temperature T e ff exhibits completely a different trend, i.e., it increases with the 
increase in centrality. The obtained values of q and T e ff for D + Au, and Au + Au collisions 
are in agreement with these tendencies. 

The contrasting behaviour of q and T e ff with respect to N part is in accord with the spirit and 
content of the non-extensive statistics, q is related with the fluctuations in temperature [Til |2T] . 
The high value of q means a high fluctuation in temperature and the system is not in the 
neighbourhood of its thermal equilibrium. The number of binary collisions in a system with 
lesser number of participant nucleons is quite low. Hence, the probability of mutual exchange 
of transverse momenta among the interacting partons in a system involving small N part is also 
quite low, which, in turn, keeps the system far away from its thermal equilibrium. On the 
other hand, an appreciable increment in the number of binary collisions is observed in a system 
possessing a large number of participant nucleons; and consequently the system can reach 
quickly to its thermal equilibrium or in the close vicinity of it. So it is quite natural that the 
fits, obtained with the non-extensive approach, will exhibit a gradual decrement in the values 
of q and increment in T e ff, as one goes from central to peripheral region or from a system with 
larger N part to that with smaller ones. 

q and T e ff change rapidly when N part lies in the range 2 < N part < 20. The change 
is rather slow for high N part . The dependence of q and T e ff on N part was observed quite 
strong over the entire N part - region while studying charged pions and proton-antiprotons[29j 
produced in nuclear collisions at ^/snn = 200 GeV. This fact is also reflected from the values 
of To and c. The values of To for charged pions and proton-antiprotons, obtained by Wilk 
et al[23] on the basis of our previous venture [29], are 220 MeV and 360 MeV respectively 
while that of c are 5.7 and 9.4 respectively. But, the values of T and c, found from the 
present analysis(T = 155 MeV and c = 1.88), are much closer to the values found from the 
analysis of pion production in e + e~ collisions [23J(T = 131MeV and c = 1.83); and T ~ 
can be recognized as the critical temperature of phase transition of the hadronizing system in 
thermal equilibrium [37]. The reason of this discrepancy in the outcome of the two studies could 
be attributed to the non- inclusion of the constraints, given by eqn(22)-(23), in our previous 
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analysis involving charged pions and proton-antiprotons. Thus, it is necessary to re-analyse 
the data on pion, proton and antiprotons produced in RHIC energies in the light of the present 
model including the constraints applied here, so that a conclusion could be reached. Lastly, but 
more importantly, there is yet another very interesting and important observation. The Fig. (2a) 
demonstrates quite convincingly that the parameter q does not seem to tend to unity for even 
higher participant numbers, so non-extensivity may be present even in the thermodynamical 
sense and not just as a simulated microcanonical effect. 

We are quite aware that so many other topical issues of interest around which we choose 
not to deal with in the present work as they demand a separate treatment. The properties of 
energy dependence of the multiplicity-density, inclusive cross-section and of the average trans- 
verse momentum behaviour are just a few among the many. They are both important and 
interesting as there are many reported anomalies and puzzles EH ED] which would be the 
subject of matter of our future studies. In fact, we commented in just the preceding paragraph 
that we would rework in the near future on the pion, kaon and proton-antiproton spectra with 
the specific concerns of constraining the parameters used by the method of calculation in the 
non-extensive approach and try to throw light on the problems and puzzles associated with the 
production characteristics of the secondaries in the high energy collisions under studies. 
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Table 1: Values of fitted parameters with respect to experimental data on ^-spectra produced 
in P + P collision at ^/s^n = 200GeV 



Np ar t 


n 


Q 


T e// (GeV) 


X 2 /ndf 


2 


0.075 ±0.003 


1.115 ±0.005 


0.097 ±0.006 


0.691/15 



Table 2: Values of fitted parameters with respect to experimental data on 77-spectra at different 
centralities of Au + Au and D ± Au collisions at RH1C 





Centrality 


Npart 


N 


c 


To(GeV) 


X 2 /ndf 




0-20 


279.9 


24.8 ±0.1 


1.81 ±0.02 


0.165 ±0.003 


0.949/7 


Au + Au 


20-60 


100.2 


9.00 ±0.03 


1.90 ±0.03 


0.160 ±0.004 


0.268/7 




60-92 


14.5 


1.31 ±0.01 


1.82 ±0.04 


0.150 ±0.005 


0.313/3 




Min. Bias 


109.1 


9.70 ±0.03 


1.83 ±0.06 


0.162 ±0.004 


0.563/7 




0-20 


15.6 


1.43 ±0.01 


1.90 ±0.02 


0.160 ±0.004 


0.580/10 




20-40 


11.1 


1.02 ±0.01 


1.92 ±0.03 


0.159 ±0.005 


0.553/10 


D + Au 


40-60 


7.7 


71 ±0.01 


1.88 ±0.03 


0.150 ±0.004 


0.237/10 




60-88 


4.2 


0.39 ±0.02 


1.91 ±0.02 


0.140 ±0.006 


1.657/9 




Min. Bias 


9.1 


0.84 ±0.04 


1.91 ±0.04 


0.151 ±0.001 


0.560/14 



Table 3: Calculated values of q and T e ff for different participant nucleons(iV part ) 



Npart 


Q 


T e// (GeV) 


2 


1.115 


0.097 


4.2 


1.101 


0.115 


7.7 


1.096 


0.123 


9.1 


1.100 


0.121 


11.1 


1.097 


0.130 


14.5 


1.093 


0.125 


15.6 


1.096 


0.131 


100.2 


1.087 


0.134 


109.1 


1.086 


0.137 


279.9 


1.085 


0.140 



14 




p T [GeV/c] p T [GeV/c] 

(c) (d) 

Figure 1: Plots of transverse momentum spectra of //-mesons produced in P + P, D + Au and 
Au + Au collisions at ^snn = 200 GeV at different centralities. The filled symbols represent 
the experimental data pointsjH]. The solid curves provide the fits on the basis of nonextensive 
approach(eqn.(21)). 
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Figure 2: Plots of the nonextensive parameter q and the effective temperature T e ff as a function 
of number of participant nucleons in Au + Au, D + Au and P + P collisions at ^s^n = 200 
GeV for production of secondary 77-mesons. 
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Figure 3: Plots of the multiplicity of r\- mesons produced per pair of participant nucleons as 
a function of number of participant nucleons in Au + Au, D + Au and P + P collisions at 
= 200 GeV. 
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